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Abstract: We present a unified treatment of exact solutions for a class of four quantum 
' mechanical models, namely the anharmonic singular potential, the generalized quantum 

isotonic oscillator, the soft-core Coulomb potential, and the non-polynomially modified 
, oscillator. We show that all four cases are reducible to the same basic ordinary differential 

Q_^| equation, which is quasi-exactly solvable. A systematic and closed form solution to the 

basic equation is obtained via the Bethe ansatz method. Using the result, general exact 
expressions for the energies and the allowed potential parameters are given explicitly for 
each of the four cases in terms of the roots of a set of algebraic equations. A hidden sl{2) 
algebraic structure is also discovered in these models. 
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1 Introduction 

Over the years, much efforts have been devoted to the determination of closed form so- 
lutions to (Schrodinger) differential equations of the form = E^, where H is the 
a I Hamiltonian and E the eigenvalue. In certain cases, the concerned equation can be re- 

duced via suitable substitutions and transformations to a well known differential equation, 
from which solutions to the original problem can be easily obtained. For instance, closed 
form solutions to many exactly solvable models in quantum mechanics have been studied 
in connection with some well known classical differential equations. However, due to lim- 
ited applications of exactly solvable systems in quantum mechanics, recent attentions have 
been on systems with partially, algebraically solvable spectra. Such systems are said to 
be quasi-exactly solvable. Thus a quantum mechanical system is called quasi-exactly solv- 
able if only a finite number of eigenvalues and corresponding eigenvectors can be obtained 
exactly through algebraic means [H El [3] . 

* d . agb oola@mat hs .uq.edu.au 
^yzz@maths. uq.edu.au 



1 



An essential feature of a quasi-exactly solvable system is that the coefficients of the 
power series solutions to the underlying differential equations satisfy three- or more step 
recursion relations, in contrast to the two-step recursions for the exactly solvable cases. 
The complexity of the three- or more step recursion relations makes it very hard (if not 
impossible) to get exact power series solutions of such systems. However, one can terminate 
the series at certain power of the variable by imposing certain constraints on the system 
parameters. By so doing, exact (polynomial) solutions to the system can be obtained, but 
only for certain energies and for special values of the parameters of the problem. 

Solutions to quasi-exactly solvable systems have mostly been discussed in terms of 
the recursion relations of the power series coefficients. However, this approach does not 
generally allow one to give explicit, closed form algebraic expressions for the allowed 
potential parameters of the systems. In this paper, we give a systematic and unified 
algebraic treatment to a class of four quantum mechanical systems, namely, 

loo e d 



1. Singular anharmonic potential [U [3 : V(r) = -uP'r^ H — j H — ^, 

2. Generalized quantum isotonic oscillator |6| : V(r) = -uP'r'^ + g-r^ 

2 (r^ -|- a'^Y 

G F 

3. Soft core Coulomb potential [7J : V{r) 



r r + j3 

^ 2 2 



loo Ar2 



4. Non polynomially modified oscillator [8] — [15] : V{r) = -co r + Y+~^~2' 

Model 1 was investigated in [3] using the Laurent series ansatz and continued fractions, 
and its ground state was found in [5]. Solutions to models 2, 3 and 4 were studied in 
jl6l-j20j on a case by case basis by means of the recursion relations. We will show that all 
four models are reducible to the same basic differential equation, which is quasi-exactly 
solvable. We solve these models exactly by using the functional Bethe ansatz method in 
|21j . Our method allows us to obtain the explicit, closed form expressions of the energies 
and the allowed potential parameters for all cases once and for all in terms of the roots of 
the algebraic (Bethe ansatz) equations. Particularly interesting are our exact results for 
model 1, in that to our knowledge this model was not previously recognized to be quasi- 
exactly solvable. We also find a underlying sl{2) algebraic structure in all these models, 
which is responsible for the quasi-exact solvability. 

The work is organized as follows. In section 2, we define the four models and transform 
the corresponding differential equations into the same basic form. In section 3, we present 
exact polynomial solutions to the basic equation underlying the four cases by using the 
Bethe ansatz method. We then give the general, closed form expressions of both the 
energies and the allowed parameters for each of the four cases. Hidden sl{2) algebraic 
structures underlying the four models are presented in section 4. In section 5 we provide 
some concluding remarks. 



2 The four models and the underlying differential equations 

In this section, we give a description of the four models to be considered in this paper. 
For each case, we reduce the underlying differential equation to the basic equation which 
is quasi-exactly solvable. 



2 



2.1 Singular anharmonic potential 

We consider a quantum mechanical model with the singular anharmonic potential [HIS], 



V{r) = \u?r^ + ^ + ^ 



(2.11 



where w, e, d > are constant parameters. The corresponding radial Schrodinger equation 
is given by 



e + l) 2 2 2e 2d 



Y{r) + 2EY{r) = 0, 



(2.2) 



where i = —1,0,1, ■■■ , and E is the energy eigenvalue. 
Making the substitution, 



Y{r) = r 



3/2+e/V2d 



exp 



UJ 2 

-— r 
2 



2d 1 
"2~^ 



yir), 



(2.3) 



we obtain 



/// X 2 / 2 3 e V2d\ , 



2ojV2d + + 



2d 



S-tj 2 + 



+ 1 



y{r) 



(2.4) 



Then the change of variable t = r'^ transforms the above differential equation into the 
form, 



-ujt'' + 2 + 



t + V2d 



y'it) + 



t 



E-uji2 + 



y{t) 



2ujV2d+ (^ + 



2d 



+ 1 



y{t). 



(2.5) 



2.2 Generalized quantum isotonic oscillator 

The generalized quantum isotonic oscillator has recently attracted a number interests. 
This system is interesting because it is endowed with properties closely related to those of 
the harmonic oscillator. The radial Schrodinger equation reads 



d^-^{r) 
dr^ 



l(i + l) 2 2 



(r2 + 0^)2 



^'(r) + 2E'${r) = 0, 



where lo, g > and a are constant parameters of the systems; £ = —1,0, 1, ■ 
the energy eigenvalue. Making the variable change z = wr^, ()2.6p becomes 



^(i2^'(z) ^ Id^(z) 



£ £(£ + 1) g(z-6ja^) 
4^ 4z ^2(z + wa2)2 



^(z) 



2 dz2 

This equation can be transformed, by the substitution. 



'2uj 



(2.6) 
and is 

(2.7) 
(2.8) 
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into the form 



z'iIj"{z) + 



7 



+ 26 + . 



2{b + l)uja^ 
z + wa^ 



i^'iz) + 



g/2 + {b + l)ii + u;a^ + 3/2) 



E 



z + uja 
7 



where 



+ 6H h 

2a; 2 4 



1 1 



(2.9) 



5 = 6(6+1) ^ 6 = 6- = -^-^^45 + 1- (2.10) 



49 + 1-1 



We find it convenient to work with the new variable t = z + ua?. Then Eq. (|2.9p becomes 
t{t - uja^)i:"{t) + -t^ + ~ \/4c/ + 1 + £ + wa^^ t + i 



-| + ^(V49 + l-l 



+ a;a" + - 



V'(i)- 
(2.11) 



2.3 Soft-core Coulomb potential 

We now consider the soft-core Coulomb potential [7] 

V{r) = ? ^ 



where G ^ Z > Q and /3 > are constant parameters. Such potential is of interest in 
atomic and molecular physics. The G = case simulates the field of a smeared charge 
and is useful in describing mesonic atoms. Solutions to the G = case has been studied in 
|19) using the method of solving the recursive relations. The parameter /3 can be related 
to the strength of a laser field, within the range of /3 = 20 - 40 covering the experimental 
laser field strength [22] . 

The corresponding reduced Schrodinger equation is 



+ 1) 2G 



2Z 



^{r) + 2E^{r) = 0, 



(2.13) 



where ii^ < is the energy eigenvalue and £ = — 1,0, 1,---. Making the substitution 

$(r) = (r + /?) r^+^e-''^''+^U{r), (2.14) 

we have 



<p"{r) + 2 



Setting 



£ + 1 ^ 1 



r + /3 



0'(r) + 2 



(£ + l)(c-l//3) + G , Z-c-{^+l)/P 



+ 



r + /3 



(2^ + c2)0(r). 

1 o 



E 



(2.15) 
(2.16) 



we obtain 

t{t + p)(t)"{t) + 2 [-ct^ + (-C/3 + ^ + 2)t + (£ + 2)/3] 0'(t) + 2t[Z -G-c{l + 2)] <P{t) 

= 2[{e+l)f3c + f3G-e-l]cl){t), (2.17) 
where we have introduced the variable t = r for later convenience. 
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2.4 Non-polynomially modified oscillator 



The non-polynomially modified oscillator potential has been found to be useful in several 
aspects of physics. In laser physics, it arises out of the Fokker-Planck equation for a single- 
mode laser |^. In field theory it provides a simple zero-dimensional model possessing a 
non-polynomial Lagrangian [9]. The 6 > case is useful in toroidal plasma [14j. 
The corresponding radial Schrodinger equation reads 



cfEjr) 



+ u; + 2A 



where a;,5 > and A and constant parameters; £ = 
With the variable change z = ujr"^, Eq. (|2.18|) becomes 



r.{r) + 2EE{r) = (2.18) 
T, 0, 1, • • • , and E is the energy. 



d'^Ejz) 1 dE{z) _ 1 
^ dz^ ^ 2~dz 4 

Making the transformation, 

E{z) 

Eq. ()2.19p becomes 

zi"{z) +\-z + 



m + l) 2A 
z + - + 



uj uj + 5z 



2uj ^ ' 







z + 



UJ 



\ i+1 


1 


j z 2 exp 


z 

2 



2z 



z + uj/6 

_ 1 
~ 2 



E X 7 



e(^) 



X/25^ +uj/6 + e + ^ 



(2.19) 



(2.20) 



z + uj/6 



(2.21) 



It is more convenient to work with the new variable t = z + uj/5. Then Eq. (|2.21|) becomes 



7\ 2uj 
+ 2 



m + 5 



E 

UJ 



uj5 



m 



A UJ 2> 
452+1+^+2 



m- (2.22) 

3 The basic equation and solutions to the four models 

It was demonstrated in the last section that the Schrodinger equations for the four cases can 
be transformed to Eqs. (|2.5|) . (|2.11|) . (|2.17p and (|2.22|) . respectively, after the appropriate 
substitutions and variable changes. These four equations have the same basic form, 



t{t - a) S"{t) + [h2t^ + hit + 5o] S'{t) + ci t Sit) = CO 5(t), 



(3.1) 



where a,b2,bi,bo,ci,co are constants. This equation has regular singular points t = 0, a 
and confluently irregular singular point t = 00. It is related, via a simple transformation, 
to the so-called generalized spheroidal wave equation (GSWE) (see e.g. [23} IM]). 

tit - to)x"it) + iBi + B2 1) x'it) + [n^it - to) - 2knit - to) + B3] xit) = 0, (3.2) 

where to,Bi,B2,B3,k,n are (possibly complex) constants. Indeed, making the substitu- 
tion Xit) = e*^*S'(t) in the generalized spheroidal wave equation yields 

t(l - to)S"it) + [2int^ + iB2 - 2into)t + Bi] S'it) + niiB2 - 2k)tSit) 
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= - {inBi + 2knto + Bs) S{t). 
This is nothing but the basic equation (|3.ip with the identifications, 



(3.3) 



to = a, n = --b2, k=-\bi + ab2- — 



Bi = bo, B2 = bi + ab2, B3 = ac\ - cq - (60 + + 0^62) • (3.4) 

Note that for the four models considered in this paper, the constants a, 621 ^i, &O7 ci, cq in 
the basic equation are real numbers. 

The basic equation (or more generally the GSWE) is quasi-exactly solvable for certain 
values of its parameters, and exact solutions are given by degree n polynomials in t with n 
being non-negative integers. In fact, the basic equation is a special case of the general 2nd 
order differential equations solved in [21] by means of the Bethe ansatz method. Applying 
the results in [21J, we have 

Proposition 3.1 Let n be any non- negative integer. The basic differential equation 
Eq. ()3.ip has degree n polynomial solutions 



5(t) = J|(t-t^), S[t) = \ for n = 



(3.5) 



i=l 



with distinct roots ti,t2, ■ ■ ■ ,tn only if the constant coefficients in (j3.1|) satisfy the con- 
straints 

ci = -nb2, (3.6) 



Co = n 



(n - 1) + 62 ^ ti + nbi, 



i=l 



where the roots ti,t2, ■ ■ ■ ,tn are determined by the Bethe ansatz equations 

2 b2tf + biti + bo . ^ ^ 

Z^T — - = — ' z = l,2,...,n. 



. tj 



ti{ti - a) 



(3.7) 



(3.8) 



In the following, we apply the above results to derive general, explicit and closed form 
expressions for the energies and the allowed potential parameters of the four models. 



3.1 Singular anharmonic potential 

In this case, S{t) = y{t), a = 0, 62 = ^1 = 2+"^^' ^0 



2d, ci 



E-uj(2 + 



and Co = T 



2u;V2d+{£ + iy-[^^ + l 



. Then by Eqs. ()3.6p and ()3.7p . we obtain 



the closed form expressions of the energies and wave functions 

e 



En = uj[2n + 2 + 



n( 

.1=1 



- ti 



exp 



UJ n 

— 
2 



2d 1 
"2~^ 



(3.9) 
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and the constraint for the potential parameters 

2 



2uj 



4n n + 1 + 



+ 



'2d J \V2d 

where the roots {tj} are determined by the Bethe ansatz equations, 



+ 1 



(3.10) 



2 + 



2d 



i = l,2,--- , 



n. 



(3.11) 



It is easy to see that y = 1 is a solution of (j2.5p provided the potential parameters 
satisfy a constraint. This solution corresponds to the n = case in the general expressions 
above. Indeed, from Eqs. (|3.9p and (|3.10p . we obtain 



Eo = uj[2 + 



Yo{r) = r3/2+-/^ 



exp 



2 



2d 1 



and the constraint 



2uV2d 



2d 



+ 1 



(3.12) 



(3.13) 



This wave function has no nodes and so the state described by it is the ground state of the 
system. Eqs. ()3.12p and ()3.13p reproduce the ground state solution and the corresponding 
constraint for the parameters found in [5J. 

For n = 1, Eqs. (|3.9p . ()3.10p and (|3.1ip give rise to the energy 



Ei=uj[A + 



and the constraint for the potential parameters 

2 



2u:[V2d + 2t^ + (^ + ^ =8 + -^ + 



2d 



+ 1 



(3.14) 



(3.15) 



where the root ti is determined by the Bethe ansatz equation, 




2d = 



2 + 



+ 4a;V2d 



It follows that the potential parameters obey 



— - 2a;V2d + 5 + 
2d 



1 

+ 2 



2 + 



+ 4a;V2d. 



The corresponding wave function is 

yi(r)=r3/2+«/^ / 2 



ti) exp 



2 



2d 1 
"2~^ 



(3.16) 



(3.17) 



(3.18) 



with ti given by ()3.16p and uj,e,d constrained by (|3.17p . This wave function gives rise to 
the first excited state of the system. 
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3.2 Generalized quantum isotonic oscillator 



In this case, S{t) = ip{t), a = ua^, 62 = -1, h = ^-^/ig'+T+i+u;a^ , bo = UJa^{^/4g'+T- 
1)' ci = Hi + - ^ - I) and Co = -g/2 + ^(V^TT -!){£ + ua'^ + 3/2). Here 

h is given in (|2.1U|) . Then by Eqs. (|3.6|) and (|3.7|) . we obtain the general, closed form 
expressions of the energies and wave functions 



En=uj{2n + l+^- ^Ag + 1^ 

n 

n 



and the constraint for the potential parameters, 



.i=l 



2,2 
r + a 



(3.19) 



n 



= -| + ^(7%+! -l)ii + u;a^ + 3/2) 
where the roots ti,t2, ■ ■ - tn obey the Bethe ansatz equations 



(3.20) 



2 _tj-{£+^-^/Af+T + uja'^)ti+uja'^{l-^/4^TT) 



. ti tj 



tiiti - uja?) 



i = 1, 2, • • • n. 



(3.21) 



It is easily seen that ■0 = 1 is a solution of ()2.11|) provided that a constraint is satisfied 
by the parameters, indeed, this solution corresponds to the n = case in the general 
expressions above. From Eqs. ()3.19p and ()3.20p . we have 



Ef) = uj 



+ ^-7%Ti) 



(3.22) 



and 



g + i + - + uja^ 



g = 2{i + l + wa^) {2£ + 3 + 2uja'^) . 
Hence the energy and wave function of the system can be written as 



and 



(3.23) 



(3.24) 



(3.25) 



The wave functions ^o(^) does not have nodes and so the state described by this is the 
ground states. 

For n = 1, Eqs. ([312]), and ^^JI^ give 



(3.26) 
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13 / 7\ 

5 + Y + 3(^ + a;a2) = 2ti + U + uja'^ + -] ^/4.g + l (3.27) 



and 



ti = - [2si + S2 + ioa^ ± x/(2si + S2)2 + uja'^i2s2 - 4si + uja'^)j , (3.28) 

where si = ^ (l - ^/4gTT) and S2 = ^ + 3/2. Substituting Eq. (^7IE\) into Eq. (fST^ . we 
have the following constraint for the allowed parameter g, 



5 + 2(£ + wa^ + 2) - + £ + wa^^ x/%TT 



+ ^ - + 1 j + c^a2 (^2£ + cjo2 + 1 + 2^4^ + 1 



which can be solved to give 



(3.29) 
(3.30) 



where 



A 



B 

C 



-151 + 93wa^ - 30£2 - QOiujo^ - 30(6^0^)^ 
-24^2^0^ - 2Ai{uja^f - - 8{uja'^f + 53 
+3 [-1380^ - 108^0^ - 1443^2 - 3246^a;a2 - 507{u}a'^f 
-2hhU^uja^ - "mUiuja^f - ^2M^ujc? - 1224{£u}a^f 

-WOSeiua^f - 6121^ - 1332(a;a2)^ - 108£^ - 300(cja^)^ 
38 + 8f + mcoa'^ + 20i + 8(uja^f + 20a;a^ 



450 



8uja^ + 8£ + 19. 
Hence we obtain the energy and wave function 

9 1 



1 B 

+ I 2A3 + — r + C 

2 3 V As 



^i(r)~/+^ (r^ + a^) 



e 2- 



r + a 



where 



+ a;a^H xiWl^H x] + ua"^ {21 + I + coa'^ + 2x) 



with 



x = 1(2AI + ^ + C 



(3.31) 



(3.32) 



(3.33) 
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3.3 Soft-core Coulomb potential 

In this case, S{t) = (p{t), a = -/3, 62 = -2c, 61 = 2(-/3c + £ + 2), bo = {i + 

ci = 2[Z-G-{£ + 2)c] and cq = 2[f3{£ + l)c + /3G - £ - 1]. Then by Eqs. and (l3TD . 

we have for the energies, wavfunction and the constraint for the parameter Z, 



Z-G 



En 



Z-G 



n 



n+£+2 2\n+£+2 
$„(r) = (r + /3) /+ie-^(''+^) 
n + 3 



i=l 

+ £-Pc\-c^ti = p{£ + l)c + pG-£-l, 



where tj satisfy the Bethe ansatz equations 



n 

Z-^ f.—f- 



cti 



-pc + £ + 2)ti-{£ + l)P 



ti{ti + l3) 



1,2,--- ,n. 



(3.34) 
(3.35) 

(3.36) 



Note that c 7^ as Z 7^ G. 

It is easily seen that </> = 1 is a solution to ()2.17p provided c and Z satisfy certain 
constraints. This solution corresponds to the n = case of our general expressions ()3.34p - 
(|3.36|) . from which we obtain the solutions 



En 



1 /Z-G 
2 



^o{r) = {r + /3)/+^ exp 

and the allowed values for the parameter Z, 

^_£ + 2 



+ 2 

£ + l-(3G 
(^ + l)/5 

G 



(r + /3) 



' + 1 



(3.37) 



(3.38) 



This wave function does not have nodes and so the state described by it is the ground 
state of the system. 

For n = 1, (I231])-(I33S1) become 



Z-G 
£ + 3 



E^ 



1 / Z-G 
'2 W + 3 



where the root ti obeys 



ch + P{£ + 2)c = 2£ + 3 - ^G, 

ctl - i-Pc + £ + 2)ti - (£ + l)/3 = 
1 



(3.39) 



(3.40) 



ti = — [-(3c + £ + 2± ^^-^c^ + 2^/3c + (£ + 2)2 
2c 



(3.41) 



It follows that the parameter Z satisfies the constraint 
(£+3)2(/3Z-^-l)(/3Z-2^-3) + (2+zy)/3(Z-G)(-3/3Z+4^+6)+2/32(Z-G)2 = 0. (3.42) 
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The corresponding wave function is 

^i{r) = (r + /3)(r - ti)/+^ exp 

where 



Z-G 

' £ + 3 



r + /3) 



(3.43) 



1 



2^(^ + 3) f{£ + 2){£ + 3) 



Z-G 



+ 



Z-G 



(3.44) 



with Z determined by p.42p . This wave function gives the first excited state of the system. 



3.4 Non-polynomially modified oscillator 



In this case, S{t) = C{t), a = f, b2 = -I, h = f + i + 2r] + I, bo = ci = 

^[§- ^-£-2r]-l] and-Co = ^ + ^+£+2rj+^. Then from Eqs. ^ and 
we obtain the energies and wave functions, 

A / 7 
En = -^ + uji2n + £+- 



En{r) ~ (1 + (5r^) /^"e^a"' 
and the allowed values for the parameter A, 



.1=1 ^ ' 



2^ 



= -(n+l) ('n + ^+£+^') 

^ ^ i=\ 



(3.45) 



(3.46) 



where the roots ti satisfy the Bethe ansatz equations 

2 _ t\-{Lolb^£^\)ti^2ujlb 

, z — i, z. 



1 7, i A 



ti{ti - u)/6) 



,n. 



(3.47) 



It is easily seen that ^ = 1 is a solution of ()2.22p provided that the parameter A obeys 
certain constraint. This solution corresponds to the n = case in the general expressions 
above. The energy and the wavefunction are 



A 



Eo = -+Lo{£ + -] , 
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(3.48) 



(3.49) 



Ho(r) ~ (1 + 5r^) r^+^e~t 
with the allowed values for A given by 

452 5 2' 

The wave functions Ho(r) does not have nodes and thus the states described by this is the 
ground state of the system. 
For n = 1, we have 

Ei = ^+uj(£ + ^Y (3.50) 
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^ h -2(^+1+'-], (3.51) 



where the root ti obeys 



452 ^ \6 2 
r, f u „ 7\ 2u ^ 



The energy and the first excited state wave function are therefore 
El = {21 + 3)(w - 45) + 45 {h " 2) - K + y j 



1 + (5r2 ti ] , (3.53) 



Ei{r) ~ (l + 5r2)/+ie-¥'^' 
where ti is given in Eq. (|3.5'2|) . 

4 Hidden Lie algebraic structure 

The basic equation (|3.1|) unifying the underlying differential equations of the four models 
considered in this paper possesses a hidden Lie sl{2) algebraic structure which is respon- 
sible for its quasi-exact solvability. 

Let us first of all recall some well-known facts. Consider the differential operators, 

J' = -r, J~^ = t^—-nt, J° = t- . 4.1 

dt' dt ' dt 2 ^ ' 

These operators satisfy the sl{2) commutation relations for any value of the parameter 
n. If n is a non- negative integer, then (j4.1|) provides a n -|- 1-dimensional irreducible 
representation, Vn+i{t) =< l,t, i^,--- > of the sl{2) algebra. From this it is evident 
that any differential operator which is a polynomial of the sl{2) generators (|4.ip will 
have the space Vn+i as its finite-dimensional invariant subspace, i.e. possesses (n + 1) 
eigenfunctions in the form of a polynomial in t of degree n. This is the main idea behind 
quasi-exact solvability of a differential operator with a sl{2) algebraization (i.e. a hidden 
sl{2) algebraic structure) [ll|2]. 

Now consider the basic equation ()3.ip and write it as the form 



HS{t) = Co S{t), H = t{t-a)^ + [h2t^ + hit + 60] ^ + ci t. (4.2) 

Then it can be checked that, if ci = —nb2 with n being any non-negative integer, the 
differential operator H above can be written as 



F = J°J°-aJ° J- +62 J+ + [n- 1 + 61] J°+ (-Ja + 60 + ? " 1 + &i • (4-3) 



2 2 



n 



2 



Thus ()4.3p provides a sl{2) algebraization for the differential operator H. In other words, 
the differential equation ()4.2p (or equivalently the the basic equation (j3.ip ) has an under- 
lying sl{2) algebraic structure □ and is quasi-exactly solvable if ci = —nb2, which is exactly 



^See e.g. [251 126) for a similar sl{2) algebraic structure in the relative motion of two charged particles 
in an external oscillator potential and the system describing charged particle moving in Coulomb and 
magnetic fields. 
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the relation (|3.6p . For such values of ci , exact solutions of the basic equation are given by 
the degree n polynomials (jS.Sp . with the allowed values of the constant coefficients 61, cq 
being determined by the constraint ()3.7p and the Bethe ansatz equations (|3.8p . 

Translating these results back to the original four cases, we can conclude that the 
four quantum mechanical models considered in this paper possess a hidden sl{2) algebraic 
structure, which is responsible for their quasi-exact solvability. 

5 Concluding remarks 

We have provided a unified treatment of the four quasi-exactly solvable quantum mechan- 
ical models and shown that the corresponding radial Schrodinger equations are reducible 
to the same basic equation which can be exactly solved by using the Bethe ansatz method. 
For each of the four cases, we have derived the closed form expressions of the energies, 
wave functions and the allowed potential parameters. We have also shown that all four 
cases possess a hidden sl{2) algebraic structure, which is responsible for the quasi-exact 
solvability of the systems. Let us remark, however, that the existence of a underlying 
Lie algebraic structure in a differential equation is only a sufficient condition for the dif- 
ferential equation to be quasi-exactly solvable. In fact there are more general (than the 
Lie-algebraically based) differential equations which do not possess a underlying Lie al- 
gebraic structure but are nevertheless quasi-exactly solvable (i.e have exact polynomial 
solutions) [21]. I 

We have seen that the basic equation ()3.ip unifying the four cases considered in this 
paper is a GSWE. Thus as by-products of the results of this paper, the GSWE possesses a 
hidden sl{2) algebraic symmetry and is quasi-exactly solvable. And its exact (polynomial) 
solutions are given by expressions ()3.5p - ()3.8p in proposition 3.1 together with the identifi- 
cations ()3.4p . It is well-known that the GSWE has important applications in astrophysics 
and molecular physics. For example, Teukolsky's equation governing perturbations of the 
Kerr black hole and the (radial or angular) equation describing the scattering of a charged 
particle on two Coulomb centers with different charges (such as one-electron diatomic 
molecules) are GSWE. It would be interesting to generalize the results of the present pa- 
per to obtain exact, closed form solutions to systems arising in the above-mentioned two 
specific physical contexts, e.g. the one-electron diatomic molecular model. Research along 
this path is underway, and results will be reported elsewhere. 
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